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Abstract 

Quantum trajectory theory is the best mathematical set up to model 
continual observations of a quantum system and feedback based on the 
observed output. Inside this framework, we study how to enhance the 
squeezing of the fluorescence light emitted by a two-level atom, stimulated 
by a coherent monochromatic laser. In the presence of a Wiseman-Milburn 
feedback scheme, based on the homodyne detection of a fraction of the 
emitted light, we analyze the squeezing dependence on the various control 
parameters. 

1 Introduction 

Photo-detection theory in continuous time has been widely developed [TJ [2J [3l 
[H [5] and applied, in particular, to the fluorescence light emitted by a two-level 
atom stimulated by a coherent monochromatic laser [SJ [5] . As well as various 
feedback schemes on the atom evolution, based on the outcoming photocurrent, 
have been proposed [7J [51 Q2] . However the introduction and the analysis of 
feedback have always been focused on the control of the atom dynamics. The 
typical aim was to drive the atom to a preassigned asymptotic state or to a 
preassigned asymptotic unitary dynamics [TU1 HH HH Q2] . 

Here, on the contrary, we arc interested in the emitted light and in employing 
control and feedback processes to enhance its squeezing properties. These can 
be checked by homodyne detection and spectral analysis of the output current. 
For these reasons we consider the mathematical description of photo-detection 
based on classical stochastic differential equations (quantum trajectories), as it 
is suitable both to consistently compute the homodyne spectrum of fluorescence 
light, and to introduce feedback and control in the mathematical formulation. 
We study how the squeezing depends on the various control parameters and 
how feedback mechanisms can be successfully introduced. We consider only 
Markovian feedback schemes a la Wiseman-Milburn [Jj [8J, as they leave the 
homodyne spectrum explicitly computable. 
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2 Detection and feedback scheme 



Consider a two-level atom with Hilbert space H — C 2 and lowering and rising 
operators er_ and cr + . Let the Pauli matrices be o~ x = er_ + c+ , a y = i(c_ — cr + ), 
= er + (7_ — er_er + and let the vector of operators (a x , o~ y ,a z ) be denoted by a. 
Let also the eigenprojectors of a z be denoted by P + = c + (T_ and P_ = (t_ct + 
and, for every angle <f>, let us introduce the unitary selfadjoint operator 

(70 = e 1 ^ er_ + c -1 ^ er + = cos </> + sin cr^ . 

A state p of the atom is represented by a point x in the Bloch sphere, 

p=^(l + x-a), fet 3 , \x\ < 1. 

We admit an open Markovian evolution for the atom, subjected to "dephas- 
ing" effects and to interactions both with a thermal bath and, via absorption 
and emission of photons, with the electromagnetic field. In particular we sup- 
pose that the atom is stimulated by a coherent monochromatic laser and that 
the emitted light is partially lost in the so called forward channel and partially 
gathered in two so called side channels for homodyne detection. 

Let the free Hamiltonian of the atom be ujo<t z /2, ujq > 0. Let the natural line- 
width of the atom be 7, let the intensities of the dephasing and thermal effects 
be given by the adimensional parameters fc<j > and n > 0, let the stimulating 
laser have frequency u > and Rabi frequency > 0. Let Auj = uo — u denote 
the detuning. 

Let the fractions of light emitted in the forward and in the two side channels 
be |ao| 2 , |ai| 2 , |a2| 2 respectively (|ao| 2 + |ai| 2 + |ck2 | 2 = 1, |ao| 2 > 0, |cni| 2 > 0, 
|«2| 2 > 0), and, for the side channels, let the initial phase of the local oscillator 
in the corresponding detector be "&k = arga^, k — 1,2. Changing $k means to 
change the measuring apparatus. Let the two homodyne photocurrents be I\ 
and Ii. 

We introduce a feedback control scheme a la Wiseman-Milburn based on 
the photocurrent I\ revealed in the side channel 1. Assuming instantaneous 
feedback, we modify the amplitude of the laser driving the atom by adding a 
term g e~ lwt I\(t) / proportional to I\, with the same frequency ui and with 
initial phase possibly different from that of the original laser. Let this difference 
be ip. 
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Then the atom has a Markovian evolution, whether we condition its state 
on continuous monitoring of the photocurrents, or we do not. Let us call a 
priori state r\ t the unconditioned one and let us call a posteriori state p t the 
conditioned one. Of course rj t is the mean of p t . Let us write the evolution 
equations in the rotating frame, where they result to be time-homogeneous. Let 
us introduce first the parameters 

c = \g\ |a |/V7 > 0, Auj c = Aiv + C7 |ai | cos(tfi - ip) e K. 

The a priori state rj t is governed by the Master equation 

drjt = Crj t At, 



Cp = — i 



+7&d {a z pa z - p)+ 7 n [a+pa^-^ {P_ , p} 



+ ~(~ 1 hil 2 ) [0--P0-+ - 2 \r. . p\ 



+7(^1 <T_-ico' ip ) p(aia + +ica v> )~ {(j«i| 2 - 2c\a\\ sin^i - <^))-P+ + c 2 , pj . 
The a posteriori state pt is governed by the non-linear stochastic Master equation 
Ap t = C Pt At + y/jV[ ai ct_ - ica v ]p t AW^t) + ^yV[a 2 a^]p t AW 2 (t), (1) 
where, for every matrix a, the superoperator T>[a] is 

T>[a]p = a p + p a* — pTi [(a + a*)p] , 

and where W\ and W2 are two independent Wiener processes. The two homo- 
dyne photocurrents are given by the generalized stochastic processes 

I k (t) = Vl\a k \Tr[o-# k p t }+W k (t). (2) 

Note that each signal term ^/7|afe| Tr [as k pt] depends on the dynamics of the 
a posteriori state pt and that, typically, it is correlated to both the noise terms 
Wi(t) and W 2 (t). 

Let us remark that, even if the feedback is based on the singular stochastic 
process I\, the mathematical formulation of the model is not affected by this 
singularity, as we do not observe directly the light in the forward channel. 

We suppose that |ao| is assigned by experimental constraints and that the 
control parameters are fi, Aw, |ai|, \a 2 \, #1, '& 2 , c and (p. Of course, if c = 0, 
then there is no feedback action on the atom, so that its a priori dynamics is 
independent of the measurement process, that is of the fractions |ai| 2 , \a 2 \ 2 and 
of the phases d 2 , ip. On the contrary, if c > 0, then the a priori dynamics is 
modified by the feedback loop and it depends also on |ai|, c and p. 
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3 Homodyne incoherent spectrum and squeez- 
ing 



We are interested in the light emitted by the atom and in particular in the 
squeezing properties of the light in the side channels 1 and 2. With the help 
of the incoherent spectrum of the homodyne photocurrents we can analyze the 
squeezing properties of the light detected in the two side channels, and thus we 
can investigate the effect of the control parameters. 

When \a-2\ 2 — 0, the fluorescence light which is not lost in the forward 
channel is gathered in a unique side channel, so that the squeezing is analyzed 
just for that light which is also detected for the feedback loop. This means that 
the eventually squeezed light would not be available for other purposes. Thus in 
this case a unique homodyne detector is employed and |ai| 2 is its efficiency. The 
meaning and the possible usefulness of the squeezing of the light involved in the 
feedback loop is discussed by Wiseman [10]. When | «2 1 2 > 0, the fluorescence 
light which is not lost in the forward channel is split in the two side channels. 
The homodyne detection of the light in channel 1 allows both the analysis of 
its squeezing and the feedback control. The light in channel 2 is detected for 
squeezing analysis, as well as it could be employed for different uses. Let us 
remark that, even if we were interested in squeezing only for channel 2, the 
choice of a feedback scheme based on homodyne detection in channel 1 would 
be still essential in order to get a time-homogeneous atomic evolution in the 
rotating frame. 

The structure of C guarantees that, for every initial preparation of the atom, 
the a priori state rj t asymptotically reaches the stationary state 



Thus we can introduce the homodyne incoherent spectrum of the light revealed 
in each side channel k as the limit of the normalized variance of the Fourier 




(3) 



where A is the 3x3 matrix giving C in the Bloch sphere language, 
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transform of the photocurrent Ik 
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It is a positive even function of its real argument p which can be computed 
from equations (p} and ((2|) by Ito calculus and by the full theory of Quantum 
Continual Measurement, which can provide the first and second moments of I\ 
[5]. Thus, for every initial state of the atom, we can obtain 



S k (vL) = l + 2 7 |a fc | 2 
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> A 2 + p 2 

where tk and s are the vectors in K 3 defined as 



tk 



(4) 



U 



Tr 



CT- Pcq 

t 2 = Tr 



e 11 p eq a+ - Tr[a i } 1 p cq ] p cq + i, — r [p eq , cr v ]) a 

\ai\ 



c M V- Poq 



Pcq e 12 <J+ - Tr[<7$ 2 p C q] Pcq) G 




More explicitly, by using the Bloch components of the equilibrium state ((3]) , we 
get 
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Each spectrum Sk depends on fed, % fl, Aoj, \ctk\, $fc, c and ip. Moreover, S 2 
depends also on |ai| and 

If Sk(p) < 1 for some /i and then the homodyne detection identified by 
"dk reveals a squeezed mode around p. of the light in channel k. 

Independently of the presence of the feedback loop, every time a parameter 
|afc| vanishes, the corresponding photocurrent Ik reduces to a pure white noise 
(shot noise due to the local oscillator) with spectrum Sk = 1 for every choice of 
the other parameters. 

Analyze first c = 0, the situation without feedback. In this case each depen- 
dence on tp disappears and S 2 becomes independent of |ai| and so that there 
is no difference between Si and S 2 . Moreover, the dependence of each spectrum 
Sk on the corresponding \ak\ reduces to the explicit multiplication coefficient 
in f4|. Therefore, when the control parameters O and Aw give squeezed light 
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in channel fc, the lowering of Sk under the shot noise level is anyhow directly 
proportional to the fraction of emitted light gathered in that channel. 

For fl = and n = there is no fluorescence light in the long run, so that 
each photocurrent Ik again reduces to a pure white noise with spectrum Sk = 1- 

For O = and n > there is no dependence on and Sk > 1. In this case 
there is only thermal light with carrier frequency ui 0l while the local oscillator 
is at frequency u>. The result are two temperature dependent Lorentzian peaks 
at /j, = ±Aw. The white noise contribution is always present. 

When Q > 0, Sk becomes ^-dependent and it can go below the shot noise 
level. This fact means that some negative correlation between signal and noise 
has been developed. Some examples are plotted for both channels, always for 
7 = 1, fed = 0, n = and |cti | 2 = \a 2 \ 2 — 0.45. Figures 1 and 2 show Si and 
5*2 respectively for Aui — (line 1) and for Aw = —2 (line 3), every time for 
values of fl and dk chosen in order to have a region with a pronounced squeezing 
(n = 0.2976, d k = -tt/2 for line 1; O = 2.0526, d k = 0.1449 for line 3). 

One could also compare the homodyne spectrum with and without n and 
fed; thus verifying that the squeezing is very sensitive to any small perturbation. 

Allow now c > 0. The optimal squeezing in channel 1 is always found for 
CI 2 = and the feedback loop is very helpful, giving good visible minima of Si 
also when |ai| 2 is not close to 1. For example, in the case \oti\ 2 — 0.45, Fig. 1 
shows Si for Auj = (line 2) and for Auj = —2 (line 4), every time for values 
of CI, c and <p chosen in order to enhance the squeezing (O = 0, c = 0.2936, 
Lp-di = ir/2 for line 2; ft = 0, c = 0.3762, i?i = 0.0482, ip = 1.9941 for line 
4). Again 7 = 1, fed = 0, n = 0. The utility of the feedback scheme can be 
appreciated by comparing lines 1 and 3 with lines 2 and 4 respectively. 

If we are interested in the light emitted in channel 2 and if |ai| 2 and \a 2 \ 2 are 
assigned by some constraints, then the squeezing in channel 2 can be enhanced 
by a feedback scheme based on the photocurrent coming from channel 1 , but the 
feedback performance is not as good as it can be for the squeezing in channel 
1 itself. Fig. 2 shows S 2 for Alu = (line 2) and for Auj = -2 (line 4), every 
time for values of ft, $1, c, ip and # 2; chosen in order to enhance the squeezing 




Figure 1: Channel 1 Figure 2: Channel 2 
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(fi = 0.2698, tfi = tt/2, c = 0.0896, 93 = 0, i9 2 = tt/ 2 f or line 2; = 2.329, 
t?i = 0.2896, c = 0.1346, p = -1.2902, ?9 2 = 0.0728 for line 4). 

Anyway, if the only constraint is |ai| 2 + |o!2 1 2 — 1 — \ a o\ 2 an d we are f ree i n 
the choice of |ai| 2 and |ct2 1 2 , then the best observable squeezing in channel 2 is 
obtained in the limit case |«i| 2 = 0, c = 0. That is, when the whole non-lost 
light is gathered just in channel 2 and the white noise I-y revealed in channel 1 
is ignored. 

Let us remark that, when we use control parameters enhancing the squeezing 
for channel k, every time Aw = we find the spectrum Sk with an absolute 
minimum in p, = 0, while whenever Aw ^ we find the spectrum Sk with two 
absolute minima, symmetric with respect to ji = 0, which turns out to be a 
local maximum. 

Finally let us remark that the idea of the papers [Til [HI H3] was to choose 
the control parameters in such a way that, in the rotating frame, the atom is 
frozen in a preassigned pure state ho G TL, i.e. in such a way that, in the rotating 
frame, both the a priori state r\ t and the a posteriori state pt asymptotically 
reach p oq = \ho)(ho\. This is possible in an exact way only in a very ideal case, 
which in our notations corresponds to Aw = 0, ip — 0, $1 = ±7r/2, |ai| = 1, 
a o = a 2 — 0, fed = 0, n — 0, which implies in particular 012 = &21 = and 
£cq = 0. Then, the a posteriori state pt is driven to a pure given state if £1 and 
c are such that y^ q + z c 2 = 1 and 2c sin $1 = 1 + z cq . But this implies t 1 = and 
the two incoherent spectra reduce to pure shot noise. This is reasonable: if the 
atom if frozen there is not incoherent scattering of light. One can check that 
actually only the coherent scattering survives, giving a ^-contribution in p = 
to the complete spectrum. If the freezing of the atom is only approximate, or if 
one tries to maximize the atomic squeezing (which is another way of stopping 
the atomic motion), one can check that all the spectra tend to become flatter 
and the squeezing tends to disappear. 
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